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SUMMARY 


A method is given for solving problems associated with Laplace and 
Poisson equations which, in general, requires considerably fewer equa- 
tions than the usual methods and which gives a convergent solution by 
the method of successive approximations. For infinite regions, by this 
method, the exact solution for the Dirichlet and Neumann problems can be 
found by solving a system of equations with as many variables as there 
are boundary points of the region. In addition, at each stage of the 
iteration a best possible estimate of the error of the approximate solu- 
tion with respect to the exact solution of the difference equation for 
the Dirichlet problem is furnished, and, for the Neumann problem, a bound 
for the error of the normal difference of the approximate solution is 
given. 


INTRODUCTION 


Certain problems in steady-state heat flow, gas dynamics, both for 
compressible and incompressible flows, plane torsion, and so forth can 
be formulated as problems associated with the Laplace or Poisson equa- 
tions in two dimensions. A frequently used method of approximating the 
solution of the Laplace equation consists of replacing the region by 
those points inside the region or on the boundary whose coordinates are 
multiples of a fixed positive number, which is the mesh size, and 
replacing the Laplace equation by the Laplace difference equation which 
says that the value of the function at a point not on the boundary is 
the mean of its values at the four neighboring points. This gives a 
system of as many equations as there are points inside the region. These 
equations are solved by relaxation or Iteration. The present report con- 
cerns a formulation of a complete system of equations for as many param- 
eters as there are boundary points of the region where the desired func- 
tion is a given linear function of these parameters for the Dirichlet 
and Ne uma nn problems. For example, for a "square region" containing, 
say, 900 inner points, the number of variables and equations required 
by the present method is 120. In any case, the larger the region the 
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greater the utility of the method. This system of equations is given 
in a form -which allows the application of the method of successive approx- 
imations, that is, each parameter is given by a linear function of all the 
parameters. In addition, the method of successive approximations applied 
to this system of equations gives a solution which is known to converge at 
least geometrically for all regions. This method is a finite-difference 
analog of the integral-equation method of potential theory and is based on 
some results of Courant, Friedrichs, and Lewy (ref. l), McCrea and Whipple 
(ref. 2), Stbhr (ref. 5)> and. Duffin (ref. 4). 

The method has the disadvantage that the calculations required for 
setting up the system of simultaneous equations are much more complex 
than in the usual methods. However, for seme problems the method can be 
used to reduce the number of variables so that the problem is within the 
capacity of an automatic computing machine. Also, if more than one prob- 
lem is to be solved for a given region this method offers a considerable 
advantage . 

Since, as is shown in this report, the Poisson difference equation 
for a finite region can readily be reduced to the Laplace equation by a 
simple computation, the above remarks apply also to the Poisson difference 
equation. 

s 

For nonlinear equations which can be written in. a form such that the 
Laplacian of the unknown function equals a given function whose arguments 
are the unknown function, its derivatives, and the space variables, the 
method of successive approximations is used, each step consisting of the 
solution of a Poisson difference equation as, for example, in the Rayleigh- 
Jansen method. 

Another advantage of the present method for the Dirichlet problem is 
that at any stage a precise bound is automatically provided for the error 
of the previous approximation and a bound is given for the corresponding 
error for the Neumann problem of the normal derivative. 

On the theoretical side the structure of the functions satisfying 
the Laplace difference equation can be completely described in terms of 
discrete potentials, and theorems which extend Duffin’ s results (ref. 4) 
can be derived. 

In addition, the method of the present report can be applied to the 
biharmonic difference equation, conformal mapping, and the theory of 
monodiffric functions (ref. 5). An interpretation in electrical network 
terms can be used to investigate electrical analog methods for solving 
the Dirichlet and Neumann problems. Finally, the results of this paper 
can be extended to n dimensions. 

b. 

In the present report the equations of the method will be derived 
and formulated for computational purposes in the analysis section. Worked , 


r, 



NACA TN 408 6 


3 


examples of the Dirichlet and Neumann methods are discussed. In appen- 
dix A the proof of Green's first , second, and third identities is given. 

In appendixes B and C the homogeneous integral-equation analogs are given 
for the Neumann and Dirichlet problems, respectively. For those readers 
who are not familiar with classical potential theory, of which the method 
presented herein is an analog, a brief sketch is given in appendix D. 

This investigation conducted at Case Institute of Technology was 
sponsored by and carried out with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


A 

A* 


(A rs ) matrix 
( A rs ') matrix 


Ar S - 9(a-rs) 

A rs' = 'P^rs ') 

Ars" = ^ark") 
a rk" = b " - c " 
( a rs) = ^yy ~ a yyt 


B ( B rs) 

B rs = 9 ( b rs) 

b 1 ( b rsO mabrix 

- (l> rs ") 

( b rs) = P 7 “ P 7 ' 


b rs' 


■u " 

b rs 


coordinates 

coordinates 


of rth point 
of rth point 


of graph of normal exterior segments 
of graph of inner normal segments 
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C 

C 

C" 


matrix analogous to kernel of integral equation of potential 
theory for Dirichlet problem 


matrix analogous to kernel of integral equation of potential 
theory for Neumann problem 


matrix used in calculation of A^W, 


"Vt 1 ^ 


c' 


D 

G 

I 


| X 7 matrix whose rows are lists of coordinates of boundary 
points 

(7 + v) x 7 matrix whose rows are lists of coordinates of 
boundary points 

diagonal matrix whose diagonal entries give number of adjacent 
inner normal points 

a region, that is, a set consisting of only inner and boundary 
points 

set of inner points of G 


1^, 7x7 identity matrix 

L = L(PjQ) = cp(P - Q) 

M 7x1 matrix each of whose rows is dipole magnitude of inner 

normal segments terminating on boundary point corresponding 
to given row 


m 7Xl matrix whose rows give mass on boundary point corre- 

sponding to row for either a simple- or double-layer 
potential 

m' 7Xl matrix whose rows give mass on normal point corre- 

sponding to row (double-layer potential) 


P 

(Pk^k) 


point or coordinates of point with integral-valued coordinates 
coordinates of kth boundary point of G 


(Pk".*k) 


coordinates of kth inner normal point of G 

coordinates of kth point of set of end points of exterior 
normal segments of G 
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') 


coordinates of kfch point of set of points consisting of 
inner normal and boundary points of G 


Q point or coordinates of point with integral-valued coordinate 


R functional operator which for given point gives value of 

function on adjacent inner points of G minus product of 
number of such points and value of function on given point 


r 


distance from origin 


s i 


ith inner normal segment 


s^ 1 ith exterior normal segment 

U function defined on G which satisfies Laplace difference 

equation on inner points 

U 7Xl matrix whose elements are values of U on boundary 

points 

V function defined on G 

W potential of simple- or double-layer distribution 

W 7Xl matrix of values of W on boundary points 

W" I x 1 matrix of values of simple-layer potential W on end 

points of exterior normal segments 

W m (v + 7) x 1 matrix of values of W on normal inner and 

boundary points 

(x,y) integral-valued coordinates 

a. 7x7 matrix where each column is list of coordinates of 

boundary points 

P number of inner normal segments 


7 number of boundary points 


7k 


diagonal elements of Dp 


factors of 7x1 matrix 
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6 difference operator 

7 Xl matrix of errors of kth approximation of U 
tj incidence matrix of graph of inner normal segments 

Tijj columns of q corresponding to normal, inner points of G 

T)p columns of q corresponding to boundary points of G 

T)’ incidence matrix of graph of exterior normal segments 

T|" columns of tj * corresponding to boundary points 

9 angle in polar coordinates 

y. magnitude of dipole for oriented unit segment 

v number of normal points 

4 number of points in graph of exterior normal segments 

a matrix whose rows correspond to normal points and which indi- 

cates boundary points adjacent to each normal point 

t number of exterior normal segments 

cp fundamental solution of Laplace difference equation 

ft constant in bound for asymptotic expression of cp 

oo Euler's constant 


Superscript : 

k denotes kth approximation of given quantity 

Subscripts: 

k denotes kth approximation of given quantity 

r,s integral values 

t transpose 

T set of boundary points of G 
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ANALYSIS 


Geometrical Definitions 

The set of points in the plane whose coordinates are integers will 
he considered. Two points will he said to he adjacent if the distance 
between them is unity. A subset of points G will he called connected 
if for ai^y two points of G there is a path consisting of segments of 
unit length connecting them in such a manner that the end points of these 
segments are all in G. A point of G is an inner point if it is adja- 
cent to four points of G. A point of G is a boundary point if it is 
not an inner point. hut is adjacent to an inner point. A point of G 
is an inner normal point if it is an inner point and is adjacent to a 
boundary point. A unit segment bounded by one inner normal point an d 
one boundary point is an inner normal segment of G. A set of points is 
a region if it is connected and every point of the set is either an inner 
or a boundary point. The symbol G will be used in the following dis- 
cussion to denote a region. The coordinates of the inner points will be 
denoted by (x,y) ; the coordinates of the inner normal points, by (p^ 1 ,^') 

where k = 1 , 2 , . . . 7 ; the coordinates of the boundary points, by 
( p k ,q -k) k = 1 , 2 , . . . v ; and the inner normal segments, by s k 

where k = 1 , 2 , . . . p. 

The matrices tj n = and T) p = (T) rik ) are defined by 


^Nik “ 



lf (v^k) 

if (P k '^ k ') 


is an end point of Sj_ 

► 

is not an end point of s-j_ 


( 1 ) 


^rik 


1 if 


is an end point of s^ 

0 if is not an end point of s^ 


( 2 ) 


It will be assumed that the boundary points of the regions considered 
consist of a finite set of points. 
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Discrete Potentials and Poisson Difference Equation 

Following the notation of SttShr (ref. 3 ) } the fundamental solution 
of the Laplace difference equation in two dimensions is written as 

<p(x,y) = q>(P) 

where x and y are integers and Cp(P) has the following properties: 

<P(x,y) = <?( I x| ,| y |) =q>(|y|,|x|) (3) 


R[cp(P)] = 



if P 
if P 


is origin *j 


is not origin 


00 


where 

R[q>(x,y)] = <p(x + l,y) +cp(x-l,y) +cp(x,y+ l) +q>(x,y- 1) -4<p(x,y) (5) 

Table I gives the values of <p(p,q)- In addition, the following asymp- 
totic estimate Is given by Stohr. There is a positive constant ft Buch 
that 





h l0 6e r 



( 6 ) 


where 


r 


= \jx? + y2 


u> = lim 
n-»<» 




In the present paper this function will be used to define a function 
of four variables. If P = (x 1 ,y 1 ) and Q = (x 2 ,y 2 ) -fcllen 

P - Q - [(x-l - x 2 ^,(y 1 - y 2 )J and L(P,Q) is defined by 


L(P,Q) = cp(P - Q) = cp (x ± - - y 2 ) 


(7) 
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By equations ( 3 ) and (4) 


L(P,Q) = cp^x-L - x 2 ),(y x - y 2 )J = cp(| - Xg| ,|yi - y 2 |) = L(Q,P) 


(8) 


and 


Hq[l(P,Q)] = Hp[L(P,Q)] = 


1 if P = Q 
0 if P ± Q 


(9) 


where Rq means that (x,y) in equation (5) is to be taken as (x 2 ,y 2 ) 
and Rp means that (x,y) in equation (5) is to be taken as (x^y-jJ. 

If n arbitrary points (x 1 ,y 1 ), (x 2 ,y 2 ), . . . (x n ,y n ^ are given 
and if to the kth point (x k ,y k ) there corresponds a real number m k 
where k = 1, 2, . . . n, then the function 


n 

W(P) = ^ mk L(P,Q k ) 


( 10 ) 


k=*l 


where Q k = ( x k ,y k ) will be called the potential of the mass distribu- 
tion whose density is m k = ni(x k ,y k ^ on the given points and zero else- 
where. By equation ( 9 ) 


R[w(P)] = 


“k lf P = Q k 
0 if P i Q k 


(H) 


where k = 1, 2, 


. n. 


If f(P) is a function defined on the inner points of a region 
then the Poisson difference equation is 


R[V(P)] = f(P) (12) 

for each inner point of the region, where V(P) may be subject to addi- 
tional conditions at boundary points and normal points. If U(P) is 
defined by 
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u(p) = v(p) - 


Y f ^k)L(p.Qk) 

k 


(13) 


where the sum is to he taken over all inner points Q^, Q 2 > • • •> then 
by equations (ll) and (12), 

r[u(p)]=o (Ik) 

at all inner points. The conditions on U(P) at the boundary points, or 
at the boundary and normal points, can be computed frcan the corresponding 
conditions on V(P) and equation (13). In this way a problem for the 
Poisson difference equation can be reduced to a problem for the Laplace 
difference equation by a direct calculation provided that the region is 
finite. For an infinite region the above reduction requires further 
study relating to the convergence of approximations for the Bum in equa- 
tion ( 13 ). 


Simple- and Double -Layer Potentials 

A simple-layer potential is the potential of a mass distribution on 
the boundary points of a region and is given by 

7 7 

W(x,y) = W(P) = £ m k L(P,Qk) = ^ m k cp [(x - p k ),(y - qj] (15) 
k^i k=l 

where m k is the mass on the kth boundary point ^ = (Pkj'lk)* Since 

there are no masses at the inner points, the simple -layer potential 
satisfies the Laplace difference equation on the inner points. It will 
be seen later that the solution of the Neumann problem can be represented 
by such a potential. 

A dipole of magnitude jj, is defined as a pair of masses situated 
on the end points of a segment where one of the masses is p and the 
other mass is -p. A double-layer potential is defined as the potential 
due to a dipole distribution on the normal segments of a region. The 
following convention will be observed: A dipole magnitude p k is asso- 

ciated with the segment s k (k = 1, 2, ... (5) and the mass of magnitude 
P k associated with s k is to be considered as being on the boundary 
point of the segment s k while the mass of magnitude -p k of this dipole 
is to be considered as being on the normal point of the segment s k . It 
Is a consequence of the definitions of and rjp (eqs. (l) and (2)) 

that 
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where k = 1, 2, 


0 

V = y %rs ^ 

lJ 

r=l 


v and 


mg 



T lrrs^ I r 


(16) 


( 17 ) 


where k = 1, 2, ... 7 and where nig and nig’ are the masses associ- 
ated with the dipole potentials at the boundary and normal points, 
respectively. 


If 

H = »1> 

is to be 


m-fc * — m~i * , mp ^ , . . . niy , — m~| , , . . . m^ , and 

Pg, . . . Pp where the subscript t means that the transpose 

taken, then the above equations can be written as 


= ri Kt p 


( 18 ) 


m = 


w 


(19) 


Hence the potential W(P) of a double layer is 


W(x,y) =w(p) = ^ m k ,( T > jj( x “ Pk* ) ,(y " ^k*)] [( x “ PkM^- Ok)] (2°) 


k=l 


k=l 


This potential satisfies the Laplace difference equation at all inner 
points which are not normal points; at normal points by equation (9), 

E [ w ( p k'’^')] "V (21) 

where k = 1 , 2 , . . . v. The double-layer potentials will be restricted 
by the condition that segments with common boundary points have the same 
dipole magnitudes . Thus , if the order of the kth boundary point is 7^ 
(i.e., the number of normal segments on this boundary point is then 

the dipole magnitudes of each of the segments is 
defined by 
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%> = Trtir ( 22 ) 

I 

It may be noted that Dp Is a diagonal matrix whose diagonal elements 

are precisely the orders of the corresponding boundary points. Thus , 
in matrix form, the above condition can be written 

H = iTpDp -1 ^ (23) 

and it will be assumed that all double-layer potentials considered ccmply 
with this condition. If both sides of equation (23) are multiplied by 
and a is defined by 

a ~ “%t^T ^ 2 ^ 

then, by equation (l 8 ) , 

m' = -aDp \i (25) 

* 

The reason for this restriction is that It guarantees the existence 
of a function U(p) defined on G so that 

* 

U(P) = W(P) (26) 

r[u(P)] =0 ( 27 ) 

for every inner point P belonging to G. To prove this assertion let 
U be defined by 

Ut = (u 1; U 2 , . . . U 7 ) (28) 

and let U(P) = W(P) if P is an inner point and U^p k ,q k ) = U k for 

boundary points. Since W(p) satisfies the Laplace difference equation 
at all nonnormal inner points it is only necessary to verify that U 
can be determined so that equation (27) is satisfied on the normal points. 

By equations ( 21 ) and (27) a necessary and sufficient condition that U 
can be determined so that equation (27) is satisfied is that the system 
of equations 

E [ w (Pk’-V) - u (Pk'<V)] - V (29) 

where k = 1 , 2 , . . . v, have a solution. Since W and U have identi- 
cal values at all Inner points, equation (29) states that at each normal 
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inner point the sum of the differences W - U at the adjacent boundary 
points is equal to the mass at the given normal point. If W is defined 
by 


w t = (^(Pl? < 3-i)^ W (p2 * 0 . 2 )* * * * ^(■^7 ,< ^7)) (50) 

then by the definition of a (eq. (24)), the system of equations (eq. ( 29 )) 
can be written 


- tJ) = m* (31) 

By equation (25) it follows that 

W - U = -Dp" 1 * (32) 

is a solution of equation (31) which defines U. This solution which 
defines U and consequently U will be called the harmonic extension 
of W on G. Hence, by equation (23) 

p. = -Tjp(W - U) (33) 

Equation ( 33 ) is the discrete analog of the discontinuity at the 
boundary of a double-layer potential of the classical theory, and, as in 
the classical theory, a double-layer potential will be used to solve the 
Dirlchlet problem. 


Reduction of Dirlchlet Problem to Integral-Equation Analog 

The Dirlchlet problem for a bounded region consists of finding a 
function U defined on a given region which satisfies the Iaplace dif- 
ference equation on the inner points of the region and which assumes 
arbitrary prescribed values on the boundary points. 

Since by equations (l 8 ), (19), and (20) W(P) is a linear function 
of p, then by_equation (30) W is a linear function of p. In equa- 
tion (33) if U is taken as given by the prescribed values of U on the 
boundary then this system of equations can be regarded as a system of 
linear equations for p. If the system of equation (33) has a solution 
for p, then m and m' can be calculated by equations (l 8 ) and ( 19 ), 
and W can be calculated by equation (20). If U(P) Is defined by the 
condition that it equals W(P) at inner points and, on the boundary 
points, coincides with the prescribed values of the Dirichlet problem, 
then U is the solution of the Dirichlet problem. This can be seen by 
the following argument. If equation (33) is multiplied by tj^ on the 

left, then equation (31) is a consequence of equation (24) . This means 
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that U is the harmonic extension of W; that is, U satisfies the 
Laplace difference equation on the inner points of G. Since by the 
construction U assumes the prescribed values, U(P) is the desired 
solution and, hence, W represents U on the inner points of G. This 
result may be formulated in the following way. For a given value of U, 
if the system of equations (33) has a solution for p, then the double- 
layer potential defined by (i represents U on the inner points of G. 

Since equation (33) is a system of (3 equations in the {3 unknowns 
|i, the existence and, in addition, the uniqueness of the solution are 
demonstrated if it can be shown that the corresponding homogeneous system 

H = -T}pW (34) 

has only the trivial solution. It will be shown in appendix C that 
equation (34) has only the trivial solution. 

The discussion above applies also to an unbounded region if the 
function U is required to be regular at infinity and to have the value 
zero at infinity. The definition of regularity given in the appendix 
corresponds to the definition of potential theory. 


Solution of Integral-Equation Analog for 
Dirichlet Problem by Iteration 

Since the systems of equations (32) and (33) are equivalent, the 
system of equations (32) will be considered. If M is defined by 



M = Dp"^m 

(35) 

then equation ( 32 ) can be written 

- 


M = -(W - U) 

(36) 

By equations (20) and (30) , 

if A = ^A rs ) and B = ( B rs ) then 



W = Am + Bm* 

(37) 

where, with r and s = 1, 

2 f ••• 7 f 


(A r s) _ 9^1 Pr “ Ps|>|^r " ^sl) 

(38) 

and with r = 1, 2, ... y 

and s = 1, 2,- . . . v. 


01 

II 

-e 

]Pr - P s ' I'Kr “ ( ls' l) 

(39) 
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Hence, by equation (25) 

W = -CM 

(4o) 

where 

C = -^Dp - Bcr) 

(41) 

Thus, equation ( 36 ) can be written 



M = U + CM 

(42) 

If M^ k ) (where k 

= 0, 1, 2, . . .) is defined by 



m(°) = U 

(43) 


M (k+1) = fj + CM^ k ) 

(44) 


Then M< k ) is the result of the kth iteration. It is known that this 
sequence converges at least geometrical!) y and an estimate of the rate of 
convergence is being investigated. 

If m( k ) and m' (4) are defined by 


m.( k ) = DpM^ 


= -oM^ 

as in equations (24) and (25), and W^ k Hx,y) is defined by 
7 V 

W^ k )(x,y) = 7 mr^^cpj^x - p r ),(y - q r )J + ^ for') k <p|(x - P r '),(y - q. r ')J 

r=l r=l 

(45) 


then the harmonic extension 


u( k ) 


of w( k ) is 


U^(x,y) = W^(x,y) (46) 

for all inner points (x,y), and u( k ) is taken as the kth approximation 
of U(x,y) on the inner points of G. 
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A best possible bound for the error of of the kth iteration 

can be given in terms of the result for the (k + l) iteration, 
and M( k ). By equation (40), 




(vr) 


and hence the values lK k ) on the boundary points of the harmonic exten- 
sion of are defined by the relation 

- U^] = + CM^ k ^ (48) 

Comparison of this equation with equation (44) yields 


U - U^ k ) = M^ k+1 ^ - (49) 

Since U - is the difference of two solutions of the Laplace dif- 

ference equation for G, U - U' ' also satisfies the Laplace difference 
equation and by the may! mum modulus principle must attain its maximum 

on the boundary. Thus the element of U - U^ k ^ which is greatest in 
absolute value is a best possible bound for the error and this can be 
computed from equation (49). Thus at each step of the iteration proc- 
ess a bound can be computed for the error of the approximation provided 
in the preceding step. 


Outline of Calculations for Dirichlet Problem 

The calculation of the solution of the Dirichlet problem may be 
divided into three stages. The first stage consists of the calculation 
of the matrix C defined by equation (4l ) ) the second stage concerns 
the solution, or the approximation of the solution, of the system of 
equation (42) j and the final stage consists of computing the values of 
the desired function on the inner points in the following way. By equa- 
tion (35)-, equation (23) can be written 

l-t = rjpM (50) 


m 1 = — oM 


and by equations (l8) and ( 19 ) 


m = DpM 


(51) 

(52) 
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Since M is computed in the second stage, these equations give m and 
m’ and by the use of equation (40), a formula for W(P) is obtained. 
Since at inner points W and U coincide, the values of U on the 
inner points are computed by this formula. The calculations of the sec- 
ond stage may be carried out in several ways, among them the iteration 
procedure defined by equations (43) and (44). It should be noted that 
any estimate of the error based upon equation (49) is equivalent to com- 
puting one step of the iteration procedure. 

The remainder of this section concerns the calculation of the 
matrix C. It may be noted that C depends only on the geometry of the 
given region and not on the values of U and, once computed, may be 
used for any Dirichlet problem for this region. 

By equation (4l) C is calculated from A, B, a, and Dp. By 
equations (23) and (24), a and Dp are calculated from qp and q^. 

as defined by equations (l) and (2) . As a preliminary step in the cal- 
culation of A and B, two auxiliary matrices, a = a rs and b = b rB 
defined, with r and s = 1 , 2 , . . . 7 , by 


a rs = (|P r - Psl’hr - *s\) (53) 

and, with r = 1 , 2 , ... 7 and s = 1 , 2 , . . . v, by 

b rs = (|Pr - Ps’ I'Kr “ 1 b‘|) (5*0 

are calculated. These two matrices are in turn calculated from three 
other matrices a yy> ^yv f ^yv' by 



where a yy is a 7 x 7 matrix, each column being a list of the coordi- 
nates of the boundary points. 
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(Pg.qg) (p 2 ,q 2 ) 



(VS) (VM 



(57) 


where $ yv Is a y x v matrix each column again being a list of the 
coordinates of the boundary points, and 



where 3^ v ' is a y x v matrix each row of which is a list of the 
coordinates of the normal points. By equations (53) and (54) 


a = cLyy - o>77t (59) 

B = P/v “ ^7v' (^O) 

and by equations (38) and (39) , 

^ rs = ^^rs) (^• L ) 

B rs = ^rs) (62) 


The procedure of stage one may be summarized as follows: For the 

given region, number the normal inner points from 1 to v, the boundary 
points from 1 to 7, and the normal segments from 1 to p. Then from 
a list of the coordinates of the normal inner and boundary points write 
the matrices a yy (eq. (56)), p 7V (eq. (57) ), and $ yv ' (eq. (58)), 

and calculate a and b. By the use of table I for 9, calculate A 
and B by equations (6l) and (62) . The next step is to write T}p 

and Tjjj by equations (l) and (2) and calculate Dp and a by equa- 
tions (22) and (24). Alternatively Dp can be written by noting that 



NACA TN 4-086 


19 


It is a diagonal matrix such that the rth diagonal element is the number 
of normal points adjacent to the rth boundary point (r = 1 , 2 , . . . 7 ), 
and a = (dik) also be written directly by the following definition 

which is equivalent to equation (24): 

f 

1 if (p k ,<i k ) is adjacent to 

°i k = ' 

0 if (Pk-f^k) no ^ ad «5 acen ' t "t° 

From A, B, a, and Dp, C is calculated by equation (4l) . As a 
partial check on the calculation of C it is proved in the appendix 
that the Siam of the elements in each row of C is zero. 

The matrices for the region indicated In figure 1 are given in 
tables II and III. 

The successive approximations of M and the error of the succes- 
sive approximations of 5 are also given In table IV. 


(63) 


Reduction of Neumann Problem to Integral-Equation Analog 

For a finite region the data of the Neumann problem are the set of 
differences for each normal segment of the value of a function at the 
inner normal end point minus Its value at the boundary end point, and the 
Neumann problem consists of determining the value of the function, subject 
to the Laplace difference equation, on the points of the region. For an 
infinite region, the additional restriction is made that the function be 
regular at infinity. As in classical potential theory it is a consequence 
of the linearity of the Laplace difference equation that the solution, if 
It exists, is determined only up to a constant. For an Infinite region, 
the solution obtained by the method of this paper is the solution whose 
value at infinity is zero and, for a finite region since a simple-layer 
potential Is used, the solution has the property that the potential repre- 
senting this function on G is also defined on the complementary region 
and has the value zero at infinity. It will also be shown in the appen- 
dix that, as In classical potential theory, a necessary condition that 
the Ne uma nn problem have a solution is that the sum of the given differ- 
ences be zero. 

The following definitions will be required. Let 


(AiU) t = (AjU l , AjU 2 , . . . AjU 7 ) 


(64) 
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where AjU k is the sum of the differences of U for the inner normal 
segments on the kth boundary point j that is, AjU k is the sum of the 

values of U on the normal points adjacent to the kth boundary point 
minus the product of the number of normal inner points adjacent to this 
boundary point and the value of the function on this boundary point. 

The segments for which at least one end point is a boundary point and 
which are not inner normal segments of G will be called exterior nor- 
mal segments of G. Let W be a function whose domain includes the end 
points of the exterior normal segments and let 

(A=w) t = (A^, A e W 2 , . . . A e W r ) (65) 

where A e W k is the sum of the values of W on those points which are 

both end points of exterior normal segments and are adjacent to the kth 
boundary point minus the product of the number of these adjacent points 
and the value of W on the kth boundary point. A formula for the cal- 
culation of equations (64) and ( 65 ) will be given in the next section. 

It may be remarked that in the notation of equation (64) the condition 
on the sum of the differences is 

7 

£ A 1 U k = 0 (66) 

k-1 

and if W is a simple-layer potential, by equations ( 15 ) and (ll). 


A i W k + A e W k " m k 


(67) 


where k = 1 , 2 , . . .7 and A^W k is calculated in the same way that 
AjU k is calculated. Unless a statement is made to the contrary all 
simple-layer potentials considered will be subject to the condition 


i 


“k “ 0 


( 68 ) 


It will be shown in the appendix that this condition insures that 

11m W(p) = 0. Since a simple-layer potential satisfies the Laplace dif- 

P-»°° 

ference equation on all inner points, it follows that W(P) is a solu- 
tion of the Laplace difference equation for the given region. Hence A-jW 
represents the differences of a solution of the Laplace difference equa- 
tion on the inner normal segments and by equation (66) 
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7 

£ AjW k = 0 (69) 

By equations (67) , (68), and (69), 

7 

£ = 0 (70) 

k=l 


For the representation of the solution of the_Neumann problem by a 
simple-layer potential it will be assumed that AqU is given. This can 
be calculated from the usual form of the data, that is, from the differ- 
ences prescribed for the inner normal segments. However, it will be 
shown in the appendix that on the inner points U is determined by AqU 

onlyj that is, if two distinct sets of differences are given so that 
AqUq = AqU 2 , then Uq_ = U 2 on the inner points up to a constant provided 

Uq and U2 satisfy the Laplace difference equation on the inner points 
of G. 


If the simple-layer potential W represents the solution of a given 
Neumann problem then 

AqW = Aqtf (71) 

and by equation (67) 

AqU + A^ = m (72) 


where mq = mq, m 2 , . . . my and m k is the coefficient of 
q^x- p k )(y “ Q-k)J in t3le representation (eq. (15)) of W(x,y). 

Conversely, if for given values of AqU, there are 7 masses 
m^, mg, . . . la? so that the simple-layer potential of these masses 
satisfies equation (72) then by equations (67) and (72) 


AqU = AqW (73) 

and hence W computed from these masses by equation (15) represents a 
solution for the given values of AqU. In addition, to show that equa- 
tion (72) has a unique solution for given values of AqU when equa- 
tion (66) is satisfied, since the number of equations coincides with the 
number of unknowns, it is only necessary to show that the corresponding 
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homogeneous system has only the trivial solution. The proof of this 
assertion will also he found in the appendix. 


Solution of Integral Equation Analog of Neumann 
Problem by Iteration 

As a preliminary calculation a formula for equation (72) will be 
derived. Assume that there are t exterior normal segments 
(s-^', 82'* . . . s T ‘) and let these segments have a fixed arbitrarily 

assigned orientation. In addition, assume that the end points of the 
exterior normal segments are denoted by (pk'^Q-k'' ) ^ ere k = 1, 2, . . . £. 

It will also be assumed that the numbering has been carried out so that 
the boundary points are the first y points of this set, that is, 

(p k j<lk) = (Pk'ScLk") with k = 1, 2, ... y. If A* = (A rs *) and m 
and TJ" are defined by 



= ( m l> m 2> * * * II V) 

W n = (w 11 w M w ,T \ 

w — y > • • • \ 

where r = 1, 2, . . . |, s = 1, 2, ... y, and = w (Pk'%^k") 

k = 1, 2, . . . then 

W" = A'm 

Let Tj 1 = ] be the incidence matrix of the graph of the 

exterior normal segments, that is. 


W 

(75) 

(76) 
where 

(7T) 


1 

if 

(iVW) 

Is terminal point of s r * 


-1 

< 

if 

(Pk'W) 

is Initial point of s r * 

' (78) 

0 

if 

(p*'W) 

is not on s r ' (r = 1, 2, . . . tj 



k = 

1, 2, . . 

• 1) 

J 



w = 
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and let tj" = (t]^" ) be the matrix consisting of the first 7 columns 
of T]', that is, 


_ 11 
%k “ 


1 if (p^lk) is terminal point of s r ' 
-1 if (p^q^) is initial point of s r * 

0 if (Pkj^k) is n °t on s r‘ ( r = 2, 

k = l, 2 , ... 7 ) 


^ (79) 


It is a consequence of 


and by equation ( 77 ) 


where 


these definitions that 
AgW = -WW" 

AgW = C'm 

C« = -^"ii’A' 


( 80 ) 

( 81 ) 

( 82 ) 


To calculate AjW, let the points of the inner normal segments 

oriented from the Inner normal points to the boundary points be denoted 
by (p^*” ) with k = 1 , 2 , . . . v + 7 where the numbering is chosen 

so that (p^ 1 " = (Pk^k) with k = 1, 2, ... 7 . Let r\ = (t^) 

be the incidence matrix of the inner normal segments where 



if (p k '”,<i k m ) 

“ (pk"’>ik"’) 
k = 1, 2, . . . 


is terminal point of s r 

is initial point of s r 

is not on s r (r = 1 , 2 , 
v + 7 ) 










(83) 


and let T]p be the first 7 columns of q. Because of the numbering 
T}p Is the matrix defined by equation (2). Also, if A" = (A^" ) is 
defined by 
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^k" = - p k )(V" - q k )] 


(84) 

where r = 1, 2 , . . .v+7 and k = 1, 2, . . .7 and 
defined by 

W ,n 

is 

ijlll _ / TJ tn TJ f »« TJ tll\ 

" VI ,w 2 * * • • "v+7 ) 


(85) 

where 



V" = w(p k '»,q k «») 


(86) 

then by equation (15) 



W*" = A"m 


(87) 

But by equation (2), equation (83), and the definition of 



AjW = -TJ rt T]W ,M 


(88) 

and hence, by equation (87), 



AjlW = C"m 


(89) 

where 



C" = -r^qA" 


(90) 


Since equation (67) holds for arbitrary values of m. 


C* + C" = i 7 (91) 

where I is the 7x7 identity matrix. Thus, C' may be calc ula ted 
by equation (82) or from equations (90) and (91). 

By equation (8l), for given values of AjU, equation (72) can be 
written 

m = AjU + C'm (92) 

The solution of this equation gives m and consequently W(x,y) which 
is U(x,y) on the given region. 
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If this system is to be solved by iteration, let 

m(°) = A^U (93) 

m ^ k+1 ) = + c*m^ k ^ (94-) 

This system is known to converge at least geometrically and an esti- 
mate of the rate of convergence is being investigated. 


Outline of Calculations for Ne uman n Problem 

The calculations by the present method may be s umma rized as follows: 
First, C 1 is calculated; second, the solution m of the system of equa- 
tions (92) is calculated or approximated; and, third, m is used to cal- 
culate W on the region by formula (15) . The procedure for the calcula- 
tion of C 1 by two methods is given below. The double calculation of 
C* is a check for the correctness of the calculation. 

In the first method A' is calculated in the following way. Let 
b' = (bj.^.*) and c' = (crk' ) defined by 



(95) 


where b* is the | X 7 matrix, each column of which Is a list of the 
coordinates of the end points of the exterior normal segments, and 




/Pi’ 9 !) (P 2 '«2) ' ' ’ (VV\ 
(Pl’»l) (V2>°-2) ■ ■ ■ (V Q -/) 

/ 

■\(P 1 ,1 1 )(p 2 ,q 2 ) • • • (P ,4,)/ 

\ f 


(96) 
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where c* is the 4 x 7 matrix, each row of which is a list of the 
coordinates of the boundary points. Let a' be defined by 


Then by equation (74) 


a* = b* - c ' 


^rk* ~ ‘P ( a rk') 


(97) 

(98) 


If t} 1 and tj" are written by definitions ( 78 ) and (79)* C 1 is 
then calculated by equation ( 82 ) . 

In the second method for the calculation of C ' , let b" = (b^ 11 ) 


b" = 



(pi "Ml"’) 

• • ■ (n-.V'A 


(V’V) 

• • * 

(P 2 '",q 2 '") 

• • • 

. . . (p 2 "',i 2 "’) \ 

• • • • • • | 

(99) 

*" ,q «") /p ««« ,q ,,, 'l 

7+v } ^7+v ) \^7+v *^7+v / 

’ * * ( P 7+v"' V+v"')/ 



where b" Is the (7 + v) X 7 matrix, each column of which is a list 
of the coordinates of the end points of the interior normal segments, 
and 


c" = 


Sl^l) W 

(W (P 2 *l 2 ) 


(P7*l 7 )\ 

( P 7'S) 


\(Pl^l) (P 2 ^2) * • * (VS)/ 


( 100 ) 


where c" is (7 + v) x 7 matrix, each row of which is a list of the 
coordinates of the boundary points. If a" = (a^" ) is defined by 


ff -uff ^fl 

l = D - C 


(101) 


Then by equation (84) 


A , " 
A rk 




( 102 ) 
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If t) and T^p are -written by equations (2) and ( 83 ) then C" is 
calculated by equations ( 90 ) and ( 91 ) , and 

C‘ = ly - C" (103) 


A partial check on the calculation of C" is that the sum of the rows 
of C" is zero. This follows from the remark that equation ( 69 ) holds 
for arbitrary values of m, if equation ( 89 ) is used and appropriate 
values of m are chosen. A check on the calculation of C* by the 
first method is that the sum of the elements of each column of C' is 
1. This is a consequence of applying the preceding remark to equa- 
tion (103). It may be noted that C* depends only on the geometry of 
the region. Thus, after numbering the inner and exterior normal segments 
and the end points of these segments, C’ is calculated, the solution of 
the system of equations ( 92 ) can be approximated by iteration (eqs. (93) 
and (94)). The desired function is given by equation ( 15 ). To illus- 
trate the method, a Neumann problem for the region of figure 2 has been 
worked and the details of the calculations are indicated in tables V 
to VII. 

The estimate of the error used is (AjU - AjU^) where 
refers to the kth approximation u(k) . If m(k) is the kth approxima- 
tion of m, and u( k ) = W^k) is computed using equation ( 15 ) where 
m 00 is the kth approximation of U, then by equation (72) 


= m^ - Agtf^ 


(104) 


But by equations (93) and (94) 


But by equation (8l) 


Aj.fi = . C.W 


. , (>t) 

AgW = C 1 m 


(105) 


( 106 ) 


Hence 


A - A -00 (k+1) (k) 

Aj_U - ApU = m - m 


( 107 ) 
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COMMENTS ON WORKED EXAMPIES 


The purpose of the worked examples is to indicate the details of 
the calculations required. It should he noted that the present method 
should not he used for problems of this size; that is, the method is most 
useful when the number of inner points is much greater than the number of 
boundary points. 

The rate of convergence for the Dirichlet problem seems to be great 
enough for practical purposes while the rate of convergence for the 
Neumann problem would seem to indicate that a modified procedure should 
be used. A tentative modification of the integral equation analog has 
increased the rate of convergence for the Neumann problem and is being 
investigated. 

The calculation of the C' matrix was checked by the calculation 
of C" . This last matrix is not given since its negative differs from 
the C* matrix by only a constant for the diagonal elements. 

Since the values of the fundamental solution were taken to four 
decimal places, the third decimal place is not exact. 

The worked example of the Dirichlet problem is for the region indi- 
cated in figure 1, and the boundary values are given In column zero of 

table III (m(°) = ft). The example of the Neumann problem Is for the 
region in figure 3, and the boundary differences are given in table VI. 
Figures 4 and 5 show the rates of convergence for these examples. 


Case Institute of Technology, 

Cleveland, Ohio, June 18, 19^6 • 
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APPENDIX A 
GREEN'S IDENTITIES 


First and Second Identities 

The proof of Green* s first and second identities for a finite set of 
points given by C our ant, Friedrichs, and Levy (ref. 1) is indicated below 
as a convenient way to explain the notation. Let Gj_ be any set of unit 
segments and let G, not necessarily a region, be the end points of these 
segments. Let U(x,y) = U(P) be a function defined on G and let 8U 
be defined for a horizontal segment as the value of U at the right end 
point minus the value at the left end point and, for a vertical segment, 
as the value at the upper end point minus the value at the lower end 
point. If V(P) is another function defined on G, consider the sum 


y V(P)r[u(P)] 

• 1 
G 

where G below the summation sign means that the Bum is to be taken over 
all points of G and R at a boundary point is to be interpreted 
according to its definition in the list of symbols. Consider the terms 
of the sum associated with the horizontal segment whose left end point 
is P and whose right end point is Q. One term arises from each end 
point and the terms may be ordered as follows: 

V(P) [U(Q) - U(P)] + V(Q)[U(P) - U(Q)] = -[v(P) - V(Q)] [u(P) - U(Q)] 


By a similar relation for the vertical segments 

y V(p)r[u(P)] = »y (&V)(8U) (Al) 

G G^ 

where the sum on the right is to be taken over all the segments of G^. 
This is Green* s first identity. 

By symmetry U and V can be interchanged on the left to get a 
similar relation, and the difference of these equations yields Green's 
second identity 

y [vr(u) - UR(V)] = 0 
G 


(A2) 
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If G is a region, U satisfies the Laplace difference equation 
on the inner points of G and V as 1; then, equation (A2) becomes 



which in the A notation is equation (66). Hence, the condition expressed 
by this equation is a necessary condition for the existence of a solution 
of the Neumann problem. 

Also, if G is a region, U satisfies the Laplace difference equa- 
tion on the inner points and V = Uj then, "equation (Al) becomes in the 
A notation. 


U t Aj_U = - \ (BU)‘ 


If U =_0 then &U = 0 for all segments of G-^. Hence, U is a con- 
stant and since U = 0, U = 0. If AjU = 0 then U is a constant. 

It is a consequence of these statements that the solution of the Dirichlet 
problem is unique if it exists and the solution of the Neumann problem is 
unique up to a constant if it exists. For finite regions these results 
are well known. However, if the regularity condition is imposed on the 
function at infinity these statements can also be asserted for Infinite 
regions. This will be proved in the next section by Bhowing that the 
above equation holds for infinite regions if U Is regular at infinity. 


Green' s Third Identity 

The proof of Green's third identity Is given as follows: Let G 

be a finite region and let V(P,Q) = L(P,Q). Then by equation (A2) 


y u(q)Rq[l(p,q)J = y l(p,q)r[u(q)] 

Q in G Q in G 

and since G is the sum of I and T (see symbol list) 

y U(Q)Rq[l(P,Q)] = W 1 + W 2 + W 5 (A3) 

L i 


Q In I 
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where 


w l(P) = £ l(p,q)r[u(q)] 


Q in I 


W 2 (P) = ^ l(p,q)r[u(q)] 


Q in r 


W 3 (p) = - ^ u(q)r q [l(p,q)] 


Q in r 


(Ah) 


(A5) 


(A6) 


Therefore, W-^(p) may be interpreted as the potential of a mass distri- 
bution on WgCP) may be interpreted as the potential of a mass dis- 
tribution on the boundary T; and W^(P) is the potential of a dipole 
distribution on the segments of N directed from the boundary to the 
normal points of dipole density -U(Q) . By equation ( 9 ) the left side 
of equation (A3) is zero if P is not in I and U(P) if P is in 
I. Thus, for any function U(P) defined on a region G 


Wi(P) + w 2 (p) + W 5 (P) 


U(p) 

< 

0 


■N 

if P is in I 

> 


if P iB not in XJ 


(A7) 


This can be stated as follows : Any function defined on G can be repre- 

sented as the sum of a potential due to a mass distribution on the inner 
points, a potential due to a mass distribution on the boundary points, 
and a dipole distribution on the normal segments. 


Let a function U(P) be defined as harmonic on a region G if 
R[u(P)J = 0 for all inner points P of G. For such a function, by 
equation (A4) , Wp = 0 and 


W 2 (P) + w 5 (p) = 



if P is in I 
if P is not in I 


(A8) 
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For finite regions with finite boundaries equation (A8) holds pro- 
vided that 


11m r% = 0 (A9) 


uniformly with respect to 0, where r and 8 are polar coordinates. 
This is proved by considering a sequence of squares with center at the 
origin whose sides approach infinity. Consider equation (A8) applied to 
the part of G in each of these squares. Sy equations (A9) and (6) it 
is seen that the parts of V ?2 and computed on the boundary of the 

square, which is slightly altered so that the region involved is har- 
monic, goes to zero and the theorem is proved. This result may be 
strengthened as indicated below. 

If r^SU is uniformly bounded and lim U = c ■uniformly with 

r-» oo 

respect to_ 6, it will be said that U is regular at infinity. If the 
function U = U - c is considered then lim U = 0 and U is regular 

r-» oo 

at infinity. Under these hypotheses it can be shown that W 2 and 

on the squares again approach zero and for U, equation (A8) holds. In 
terms of equation (A8) this becomes - _ 


I 


l(p,q)r[u(q)| 


Q In r 


- £ u(q)r q [l(p,q)] + 

Q in r 


£ cR Q [L(P,Q)] 
Q in T 


U(P) - 


0 


c 


if P is in I 
if P is not in I 


(A10) 


To evaluate the last term on the left consider a sufficiently large 
square modified as above containing F. Let the part of G contained 
in the square be denoted by G' and let the boundary of G* be denoted 
by T + T' . If equation (A8) is applied setting U = c, 


W C + W C * 


’ c if P 
0 if P 


is in Interior of G' 
is not in interior of G* 
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where 


W C = - 


£ cHqL(p;q) 


Q in r 


¥ 1 = - 
C 


£ cRqL(P,Q) 

Q in r 

But by equation (A8) applied to the square and its interior 


v - 


= < 


c if P is in interior of square 
0 if P is not in interior of square 


Hence, if P is in both the int erica? of G* and the interior of the 
square, Wq = 0. Since the square may be taken as large enough to 
include any finite point of the interior of G, Wq' = 0. 

Hence, in equation (A10) , the last term on the left side is zero 
and equation (a 8) for harmonic functions regular at infinity becomes 


W 2 (P) + W 3 (P) + c = j 


fu(P) if P is in I 


if P is not in I 


(All) 


This formula can be generalized to the case where U is not reg ular 
at infinity. From the discussion above. 


Z 1 if P is in I 
Rn[L(P,Q)] = - 

0 if P is not in I 
Q in T ^ 

This suggests, as in potential theory, defining the mass of U as 


(A12) 


£ r[u(q2 


(A13) 


Q in T 


M(U) = - 
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If U can be represented in the form 

U(P) = ML(P,P' ) + U^P) 

■where U^P) is regular at infinity then equation (All) can be shown to 

hold. It may be noted that definition (A13) resembles the definition of 

flux in potential theory. 

If the components of U are identical then U is a constant func- 
tion and = Wg = 0 for all values of P. Also, Wj(P) can be inter- 
preted as the potential due to a dipole distribution for which all the 

components of p equal the same constant, say c. Then by equation (a8) 

such a potential is c on all inner points and zero on all other points. 
Thus, by equation (40) the sum of the columns of C must be zero. 

A similar analysis of the Green's identities for infinite regions 
applied to simple-layer potentials shows that condition (68) guarantees 
the validity of Green's identities for functions represented by simple- 
layer potentials. 
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APPENDIX B 

HOMOGENEOUS INTEGRAL-EQUATION ANALOG FOR NEUMANN PROBIPM 


The homogeneous system corresponding to equation ( 72 ) is obtained 
by setting AjU =0 in equation ( 72 ); that is. 


AgW = m 


(Bl) 


By equation ( 68 ) and the discussion of appendix A W is reg ular at 
infinity and its value at infinity is zero. Thus, for finite and infi- 
nite values of G, Green* s identity gives 


y (sw ) 2 = j 


A_. 


•w 


(B2) 


But by equations ( 67 ) and (Bl) 


AjW = 0 


Hence, by equation (B2) SW = 0 and W is a constant on G. Now on 
the set consisting of boundary points of G or the points of the com- 
plement of G, W is harmonic. Since the set of boundary points of 
the above set is also the set of boundary points r of G and since 
W is constant on P, it follows that W is constant on the complement 
of G. But the point at infinity is in one of these regions an d w is 
zero at this point. Hence, W is identically zero. Hence, AgW = 0 
and m = 0. Thus the homogeneous integral analog system for the Neumann 
problem has only the trivial solution and the existence of a unique solu- 
tion for the system of equation ( 72 ) is proved. 
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APPENDIX C 

HOMOGENEOUS INTEGRAL-EQUATION ANALOG FOR DIRICHLET PROBLEM 


The homogeneous integral-equation analog for the Dirichlet problem 
is given as follows: In order to prove the existence and uniqueness of 

the system of equations (33) or ( 36 ) it suffices to prove that the system 
of equations 


M = -W (Cl) 

obtained by setting U = 0 has only the trivial solution. 

As a preliminary step a representation of the harmonic extension 
of W will be derived. Let V be that solution of the Laplace dif- 
ference equation on G for which the boundary values V of V are 

V = M (C2) 


By equation (A8) 

Y L(P,Q)R[v(Q)] - £ V(Q)Rq[l(P,Q)] = 

Q in T Q in T 

It is a consequence of equations (C2) and (20) that 


V(P) for P in I 
0 for P In T 


(C3) 


w(p) = - £ v(q)r q [l(p,q)] 


Q In r 


If H(P) is defined by 


H(P) = £ L(P,Q)r[v(q)] 


Q in T 

then H(P) is harmonic in G, and 

r 

V(P) - H(P) for P in I 
[V(P) - H(P)] - V(P) for P in T 


W(P) = J 


(c4) 


( 05 ) 


(c6) 
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Also, since 


I 


r[v(q)] = 0 , 


Q in r 

zero at infinity. 


H(P) 


is regular at infinity and is 


From equations 


( 36 ) and (C2) 


W(P) 


U(p) for P in I 
U(P) - V(P) for P in T 


Comparison with equation 


(C 6 ) yields 
U(P) = V(P) - H(P) 


for all points 
If U = 0 


P in G. 

then by equation (c8) 


V(P) = £ m(Q b .)L(P,Q !£ ) 

9k inr 

where 

"•(%) - ®[v(%)j - AiV k 

But 

m(Q k ) = AiV k + A e V k 

Hence 

A j = 0 

e 


(C7) 


(C 8 ) 


(C9) 


(CIO) 


(oil) 


By reasoning similar to that of the- preceding section V is identically 
zero and hence M is zero. This completes the proof of the existence 
and uniqueness of the solution of the integral-equation analog for the 
Dirichlet problem. 
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APPENDIX D 

CLASSICAL POTENTIAL THEORY 

Since the theory in this paper is an analog of classical potential 
theory, a brief sketch of this subject and the reduction of the Neumann 
and Dirichlet problems to integral equations is given here. All regions 
mentioned below are assumed to be at least as regular as is required for 
the following statements to hold. 

If r is the distance from a given point to a unit positive elec- 
trical point charge, the value of the potential of the charge at the 
given point is l/r in suitable units. This potential function is a 
solution of the Laplace partial-differential equation 

a 2 f a% a% 

8x 2 ay 2 8z 2 

in any region not containing the charge. If S is a surface and if a 
charge density m 1 (x,y,z) is defined on the surface, then the poten- 
tial H-]_ of this charge distribution is 



The potential is also a solution of the Laplace differential 

equation in any regular region not containing S and is called a simple- 
layer potential. This can be verified by direct substitution if the 
derivatives are calculated by inverting the order of differentiation and 
integration. It is a central theorem of classical potential theory that 
any solution of the Laplace equation in a region can be represented as a 
simple-layer potential. 

If two numerically equal charges of opposite sign are brought together 
along a fixed line and the magnitudes of the charges are varied so that the 
product of the numerical value of the charges and the distance between them 
Is held fixed, then the limit of the potential obtained by letting the dis- 
tance between the charges go to zero is called the potential of a dipole 
oriented along the given line and has the form 8(l/r)/8n, where 8/8n 
denotes differentiation along the given line. This element is called a 
dipole and is also a solution of the Laplace differential equation. As 
above, if a dipole density mg is defined on S, the potential 
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Hg(x,y,z) = JJ |^(J)dA 


where 

r = /(x - |) 2 + (y - tj ) 2 + (z - g ) 2 

and 3/6n denotes differentiation along the normal to S, is also a 
solution of Laplace's differential equation in any regular region not 
containing S. It can he shown that any solution of the Laplace equa- 
tion can he represented, up to an additive constant, hy H 2 , that is, 
hy a dipole or double-layer potential. 

The Diriehlet or first boundary- value problem of potential theory 
consists of finding a solution of Laplace's equation in a region which 
assumes given values on the boundary of the region. The second problem 
is to find a solution in the region in which the normal derivative on 
the boundary is a given function. The first problem can be reduced to 
the problem of determining a dipole density on the boundary whose poten- 
tial is the desired function; the second problem can be formulated as 
the problem of determining a charge density on the boundary whose poten- 
tial is the desired function. In either case an integral equation for 
the dipole or charge density can. be derived by using either the discon- 
tinuity of a dipole distribution on the boundary or the discontinuity of 
the normal derivative of a simple-layer potential on the boundary. 

For the Laplace equation in two dimensions the physical interpreta- 
tion is not so straightforward, but by analogy with the three-dimensional 

case a "charge" is considered whose potential Is log e = =4 where r is 

the distance between the charge and the point at which the potential is 
being evaluated. Dipoles are as defined previously, but Instead of con- 
sidering distributions on surfaces, distributions on curves are consid- 
ered. All the other statements for the three-dimensional case hold for 
the two-dimensional case. Thus the potentials of a charge and the dipole 
distributions, respectively, are 

H 1 (x,y) ,q)log e <Is 


H 2 U,y) = J “2(6^)^ ^Log e i)ds 
^ c 
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where r = l/(x - |)^ + (y - q)^, d/dn indicates differentiation along 
the normal to the curve C, and ds is the element of arc length on C. 

If the limit from the interior of a region is denoted by the subscript i 
and the limit from the exterior by the subscript e, the discontinuity of 
a double-layer potential on the boundary of a given region can be expressed 
by 


H 2i (x,y) = H 2 (x,y) + itm^x^) 

H 2e (x,y) = H 2 (x,y) - xcm 2 (x,y) 

where (x,y) is a point on the boundary. For a simple-layer potential if, 
at a fixed boundary point (xQ,y 0 ), differentiation in the direction of 

5/cta then 

-jtm^Xojyo) + J m;L (i,Ti)|-(log e j)ds 


rtm l(Wo) + f m l^^^( log e f)* 8 

of the given region. 

In order to solve the Dlrichlet problem by representing the solution 
as a double-layer potential, it is noted that H 2i (x,y) must coincide with 
the given boundary value of the desired function. Substituting for H 2 
its representation in terms of m 2 yields 

m 2 (|, i i)|-(l o g e i^ds + = F(s) 

where F(s) is the given value of the desired function at the point on 
the boundary curve C, x = x(s), and y =~y(s}. If the substitutions 
are carried out, a linear integral equation for m 2 is obtained. The 

exterior problem leads to a similar equation if the second discontinuity 
condition for a double-layer potential is used. 


this normal is denoted by 

i 


e 

where C Is the boundary 


- SHj/x^yo) 
eta 

~aHi(x 0 ,y 0 ) 

§n 


For the inner Neumann problem it is noted that 


dn 


is given 


If the solution is regarded as represented by a simple-layer potential 
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and , hence , the first discontinuity condition for the normal derivative 
of a simple-layer potential becomes an integral equation for the "charge 
density." Sufficient conditions which ensure the existence and uniqueness 
of the solutions of these problems as well as the proofs of the above 
statements can be found in reference 6. The outer problem can be solved 
by using the second condition. 

The potential theoretical approach to the Neumann and Dirichlet 
problems may be summarized in the following way. The given problem con- 
sists of finding a function which satisfies the Laplace equation on a 
region and, in the case of the Dirichlet problem, assumes given values 
on the boundary of the region or, in the case of the Neumann problem, 
has a normal derivative on the boundary which assumes prescribed values 
on the boundary. In either case the unknown function is regarded as the 
potential of a simple- or double-layer charge density on the boundary 
which satisfies a certain integral equation. Thus the original problem 
is reduced to the problem of solving an integral equation for the charge 
density. Once the charge density is known, the potential, which is the 
desired function, can be computed directly. 

The methods given in this paper for handling the corresponding prob- 
lems for the Laplace difference equation are analogs of the above methods. 
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♦ I 


TABIf X 

IAPXACE SOIDUCH, q>(p,i) 


ES 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

1 

s 

3 

4 
3 
6 

7 

8 
9 

000,000 

250,000 

318,330 

363,380 

586,620 

424,415 

450,281 

440,376 

462,207 

488,075 

476,994 

482,595 

495,962 

513,944 

533,548 

512,902 

516,250 

525,303 

538,190 

553,152 

568,916 

342,116 

544,400 

550,811 

560,559 

571,936 

584,679 

597,853 

566,760 

568,422 

573,181 

500,480 

589,633 

599,992 

611,027 

622,358 

588,082 

589,347 

593,013 

596,746 

606,103 

614,631 

623,923 

633,630 

643,559 

606,874 

607,870 

610.778 
615,387 
621,403 

628,505 

636,386 

644.779 
653,468 
662,263 

S 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

■n 

H 

625,676 

624,480 

626,842 

650,622 

655,6X8 

641,599 

648,555 

655,613 

665,251 

671,096 

679,056 

638,869 

639,532 

641,488 

644,640 

648,847 

653,939 

659,742 

666,088 

672,824 

679,822 

686,974 

694,194 

652,735 

653,292 

654,938 

657,605 

661,190 

665,569 

670,607 

676,171 

682,136 

688,393 

694,845 

701,413 

706,053 

665,488 

665,962 

667,366 

669,651 

672,740 

676,540 

680,946 

685,852 

691,157 

696,767 

702,599 

708,581 

714,654 

720,766 

677,294 

677,702 

678,914 

680,892 

SgS 

690,784 

695,135 

699,873 

704,920 

710,203 

715,659 

721,233 

726,878 

732,555 

688,283 

688,659 

689,695 

691,424 

693,782 

696,713 

700,152 

704,031 

708,281 

712,836 

717,654 

722,620 

727,742 

732,959 

738,232 

743,331 

698,562 

698,874 

699,803 

701,327 

703,412 

706,013 

709,079 

712,555 

716,384 

720,509 

724,879 

729,443 

734,157 

738,982 

743,883 

748,831 

753,799 

706,217 

708,493 

709,317 

710,669 

712,525 

714,846 

717,597 

720,726 

724,189 

727,938 

731,926 

736,116 

740,461 

744,929 

749,488 

754,109 

758>768 

763,445 

717,319 

717,566 

718,300 

719,509 

721,171 

723,257 

725,734 

728,565 

731,709 

735,127 

738,780 

742,630 

746,643 

750,786 
TO, 030 
759,348 
763,719 
768,122 

772,540 

725,926 

726,150 

726,809 

727,896 

729,593 

731,276 

733,519 

736,089 

738,954 

742,080 

745,434 

748,983 

752,695 

756,542 

760,497 

764,556 

768,638 

772,784 

776,957 

781,143 


4 =" 

VjJ 
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TABIE I.- Continued 
LAPUCE SOLUTION, 9<p,q) 


V 

4 

20 

21 

22 

23 

2 k 

£5 

26 

27 

28 

0 

734,096 

741,864 

749,271 

756,348 

763,123 

769,622 

775,866 

781,874 

787,663 

1 

734,295 

742,045 

749,435 

756,496 

763,362 

769,750 

775,964 

781,983 

787,7® 

2 

734,890 

742,585 

743,477 

749,928 

756,949 

763,675 

770,131 

776,336 

782,310 

788,069 

3 

735,873 

750,741 

757,694 

764,360 

770,763 

776,921 

782,652 

788,373 

4 

757,233 

744,709 

751,866 

758,725 

765,309 

771,638 

777,731 

» 

789,274 

3 

738,936 

746,265 

753,290 

760,032 

766,312 

772,750 

778,762 

790,166 

6 

740,975 

748,127 

754,995 

761,600 

767,959 

774,088 

780,003 

783,717 

791,242 

7 

743,318 

750,271 

756,964 

763,414 

769,634 

775,641 

781,443 

787,060 

792,496 

8 

745,939 

752,675 

739,177 

765,456 

771,325 

777,395 

783,077 

788,582 

793,918 

9 

748,807 

755,314 

761,612 

767,709 

773,614 

779,338 

784,887 

790,272 

795,499 

10 

751,894 

758,163 

764,247 

770,153 

775,886 

781,433 

786,862 

792,U9 

797,230 

n 

755,171 

761,197 

767,061 

772,769 

778,323 

783,728 

788,98 9 

794,112 

799,301 

12 

758,610. 

764,391 

770,033 

775,538 

780,909 

786,147 

791,255 

796,238 

801,100 

13 

762,189 

767,723 

773,142 

778,443 

783,627 

788,695 

793,647 

798,487 

803,217 

Uv 

765,880 

771,171 

776,368 

781,465 

786,462 

791,357 

796,152 

800,846 

803,441 

15 

769,661 

774,715 

779,692 

784,587 

799,596 

792,^20 

794,121 

798,737 

803,304 

807,763 

16 

773,514 

778,336 

783,098 

787,794 

796,973 

801,450 

803,830 

810,172 

17 

777,420 

782,018 

786,571 

791,072 

795,516 

799,900 

804,219 

808, 4p 
811,164 

812,65? 

815,214 

18 

781,362 

785,744 

790,095 

794,406 

798,673 

802,891 

807,055 

19 

785,328 

789,503 

793,658 

797,785 

801,880 

800,935 

809,947 

813,912 

817,82? 

20 

789,304 

793,281 

797,248 

801,198 

805,125 

809,022 

812,883 

816,710 

820,494 

£1 

22 

23 

24 

as 

28 

27 

28 

797,068 

800,835 

804,471 

804,635 

808,086 

811,544 

808,399 

811,694 

815,003 

B18,317 

812,143 

815,290 

818,435 

821,631 

824,813 

815,661 

818,867 

821,895 

824,940 

827,993 

831,054 

819,549 
822,421 
823,320 
828,239 
831,172 
634, 114 
837,060 

823,204 

825,949 

828,723 

831,525 

834,340 

842,848 

X 

29 

30 

31 

32 

33 

34 

35 

36 

3T 

0 

793,249 

793,344 

796,646 

803,866 

808,919 

813,818 

818,370 

823,184 

827,668 

832,029 

1 

756,733 

803,949 

806,997 

815,891 

818,639 

823,24? 

823,444 

827,729 

832,087 

2 

793,628 

796,999 

804,197 

809,230 

814,110 

818,843 

827,913 

832,262 

3 

794,098 

799,439 

804,609 

809,617 

814,474 

819,188 

823,767 

828,219 

832,551 

4 

794,752 

800,051 

805,182 

810,155 

814,980 

819,663 

824,218 

628,643 

832,953 

833,471 

5 

795,384 

800,830 

805,912 

810,842 

815,626 

820,274 

824,793 

829,190 

6 

796,590 

801,772 

806,796 

811,673 

816,409 

821,013 

825,490- 

829,®! 

834,097 

7 

797,763 

802,871 

807,829. 

812,644 

817,324 

821,877 

826,308 

830,624 

834,830 

8 

799,095 

804,121 

809,003 

813,750 

818,368 

822,863 

827,241 

851,508 

833,669 

9 

800,577 

805,513 

810,314 

814,9® 

819,534 

823,965 

828,284 

852,497 

836,606 

10 

802,202 

807,041 

811,733 

816,345 

820,817 

825,179 

829,435 

835,589 

857,645 

11 

803,961 

808,697 

813,314 

817,817 

822,211 

826,500 

830,687 

S8:K 

838,773 

12 

805,642 

810,471 

814,989 

819,401 

823,711 

827,921 

832,036 

839,995 

13 

807,838 

812,356 

836,771 

821,088 

825,309 

829,438 

833,478 

837,431 

841,300 

1*- 

809,939 

814,342 

818,651 

822,870 

827,000 

831,044 

835,005 

836,614 

638,885 

842,686 

15 

612,133 

816,421 

820,623 

824,741 

826,777 

832,734 

840,417 

844,148 

16 

8l4,4lfl 

818,385 

822,677 

826,693 

830,634 

834,502 

838,2S6 

842,024 

845,681 

17 

816,777 

820,826 

824,807 

826,719 

832,563 

836,341 

840,052 

843,699 

847,282 

18 

819,205 

823,136 

827,005 

830,813 

834,560 

838,246 

841,872 

845,438 

848,943 

19 

821,694 

825,506 

829,265 

832,969 

836,618 

840,212 

843, pi 
845,685 

847,235 

850,666 

20 

824,255 

827,931 

830,402 

831,579 

835,179 

838,730 

842,232 

849,087 

852,440 

21 

826,822 

833,941 

837,438 

840,892 

844,302 

847,668 

830,988 

®4,264 

22 

829,448 

832,914 

836,343 

839,741 

843,098 

846,417 

849,696 

852,935 

856,153 

858,043 

23 

832,106 

635,460 

838,786 

842,001 

845,543 

848,572 

831,765 

854,922 

24 

834,790 

838,036 

841,258 

844,454 

847,622 

850,761 

853,869 

®6,946 

®9, 990 

25 

837,495 

840,216 

840,634 

843,755 

846,654 

849,930 

832,981 

856,006 

859,002 

861,970 

26 

843,252 

846,274 

849,278 

852,264 

855,228 

858,170 

861,087 

863,980 

27 

842,949 

843,884 

848,809 

851,721 

854,618 

857,496 

860,358 

863,196 

866,016 

28 

845,689 

848,526 

851,358 

654,180 

856,990 

859,786 

862,567 

8®, 931 
867,482 

868,075 

29 

848,432 

851,175 

833, 9!5 
856,479 
859,045 

656,650 

859,376 

862,091 

864,794 

870,155 

30 

31 

32 

33 
3* 

35 

36 

37 

853,827 

859,126 

861,612 

864,098 

861,772 

864,175 

866,584 

868,995 

864,408 

866,734 

869,068 

871,406 

873,746 

867,034 
869,286 
871, 54T 
873,815 
876,086 
878,359 

869,649 

871,850 

674,020 

876,21? 

878,424 

880,632 

882,842 

872,251 

SS:iS 

878,61a 
880,758 
682,903 
6®, 052 
887, 003 
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TABLE II.- MATRICES FOR DIRICHLET PROBLEM - Continued 

(c) a matrix 


Column 


2 5 


8 9 10 11 12 13 14 


1 0,0 1,0 
2 1,0 0,0 
3 2,0 1,0 

4 3,0 2,0 

5 4,1 3,1 

6 4,2 3,2 

7 4,3 3,3 

8 3,4 2,4 

9 2,4 1,4 

10 1,4 0,4 

11 0,4 1,4 

12 1,3 2,3 

13 1,1 2,2 

14 1,1 2,1 


2.4 1,4 

1.4 0,4 
0,4 1,4 

1,4 2,4 
2,3 3,3 
2,2 3,2 
2,1 3,1 
1,0 2,0 

1,0 


1.3 1,2 

2.3 2,2 

3.3 3,2 

4.3 4,2 
5,2 5,1 

5,1 5,0 

5.0 5,1 

4.1 4,2 

3.1 3,2 

2.1 2,2 
1,2 
0,1 
0,0 


(d) b matrix 


Column 


2 

3 

1,1 

2,1 

0,1 

1,1 

1,1 

0,1 

2,1 

1,1 

3,0 

2,0 

3,1 

2,1 

3,2 

2,2 

2,3 

1,3 

1,3 

0,3 

0,3 

1,3 

1,3 

2,3 

2,2 

3,2 

2,1 

3,1 

2,0 

3,0 


9 10 


0,2 












































TABI£ II.- MA 35 UCBS POft DIRICHIET PROBiaJ - Continued 
(o) A matrix 


Column 


2 

> 

mm 

5 

6 

0.2500 

0.3634 

0.4303 

0.4824 

0.4960 

0 

.2500 

.3634 

.44o4 

.4622 

.2500 

0 

.2500 

.3866 

.4244 

.2634 

.2500 

0 

■3183 

.3866 

.4404 

.3866 

.3185 

0 

.2500 

.4622 

.4244 

.3866 

.2300 

0 

.mi 

.4622 

.4404 

.3634 

.2300 

.4960 

.4824 

.4770 

.44o4 

.3866 

.4824 

■ 4770 

.4824 

.4622 

.4244 

■4770 

.4824 

A960 

.4881 

.4622 

.4824 

.4960 

.5140 

■3l4o 

.4960 

.4622 

.4881 

.5140 

.5253 

.5163 

.4244 

.4622 

.4960 

.5163 

.5129 

.3866 

.44o4 

.4824 

.5129 

.5163 


0 . 53^0 

.4881 

.4822 

.4404 

.5654 

.2500 

0 

.3183 

xfljC £ 

,pooo 

.4404 

.48a4 

.5129 

.3163 

.5233 


8 

9 

10 

11 

-1 

12 

0.5140 

0.4960 

0.4824 

0.4770 

0.44o4 

.4960 

.4824 

.4770 

.4824 

.4622 

.4824 

.4770 

.4824 

.4960 

.4801 

.4770 

.4824 

.4960 

.5140 

-5140 

.4404 

.4622 

.4881 

.5140 

■5£53 

.3866 

.4244 

.4622 

.4960 

.5163 

.3183 

.3866 

.4404 

.4824 

.512? 

0 

.2300 

.3634 

.4303 

.4824 

.2500 

0 

.2500 

.3634 

.44o4 

.3634 

.2500 

0 

.2300 

.3866 

.4303 

.3634 

.2500 

0 

.3183 

.4824 

,44o4 

.3866 

.3183 

0 

.4960 

.4622 

.4244 

.3866 

.2500 

.5140 

.4881 

.4622 

•44o4 

.3634 


0.3866 

.4244 

.4622 

.4960 

.5163 

■5129 

.5163 

.4960 

.4622 

.4244 

.3866 

.2300 

0 

.2300 


0.3183 

.3866 

.44o4 

.4824 

■5129 

.5163 

.5253 

.5140 

.4881 

.4622 1 

.4404 

-3654 

.2300 

0 


(f) B matrix 


Column 


H 

2 

3 

-r 1 

5 

6 

7 

8 

9 

10 

0.2500 

0.3183 

0.3866 

0.4404 

0.4622 

0.4881 

0.4622 

0.44o4 

0.4503 

0.3634 

•3183 

.2500 

•3183 

.3866 

.4244 

-.4622 

.4404 

.4303 

.4404 

.3866 

.3866 

.3183 

.2500 

.3183 

.3866 

.4404 

.4303 

,4404 

.4622 

.4244 

.4404 

.3866 

.3183 

.2500 

.3654 

.4303 

.4404 

.4622 

.4881 

.4622 

.4770 

.4303 

.3634 

.2300 

*3183 

.3866 

.4244 

.4622 

.4960 

.4824 

.4824 

.4404 

.3866 

-3183 

.2500 

.3183 

.3866 

.4404 

.4824 

.4770 

.4960 

.4622 

.4244 

.3866 

.3183 

.2500 

-3634 

.4303 

.4770 

.4824 

.4881 

.462? 

.44o4 

.4303 

.3634 

.2500 

-3183 

.3 866 

.4404 

.4622 

.4622 

.4404 


.4404 

.3866 

.3183 

.2500 

.3183 

.3866 

.4244 

.4404 

.4303 


.4622 

.4244 

.3866 

.3183 

.2500 

-5183 

.3866 

.4303 

.4404 

.4622 

.4881 

.4622 

.44o4 

.3866 

•5185 

.2500 

.3634 

.3866 

.4244 

.4622 

.4960 

.4824 

•4770 

.4303 

*3634 

.2500 

.3183 

.3183 

.3866 

.44o4 

.4824 

.4770 

.4824 

.4404 

.3866 

.3183 

.2500 

.2500 

.3634 

•4303 

.4770 

.4824 

.4960 

.4622 

.4244 

.3866 

■3183 
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TABLE II.- MATRICES FOR DIRICHIET PROBIBrf - Concluded 


(g) C matrix 


Row 

Column 

1 

JL 

r\ 

C. 

n . 

J 

1, 

T 

c 

C 

u 

•■7 

1 

G 

u 

rt 

7 

1 A 
-LW 

“i 1 

JUL 

*i r\ 
JL C. 

1 1 

14 

1 

0.2500 

0.0683 

0.0232 

0.0101 

- 0.0420 

-0.0338 

- 0.0259 

-0.0259 

-0.0338 

- 0.0420 

- 0.0467 

-0.0101 

-0.0232 

- 0.0683 

2 

.0685 

.2500 

.0683 

.0232 

-.0538 

-.0378 

-.0259 

-.0338 

-.0420 

-.0467 

-.0420 

-.0233 

-.0378 

-.0683 

5 

.0252 

.0683 

.2500 

.0683 

-.0683 

-.0378 

-.0218 

-.0420 

-.0467 

-.0420 

-.0338 

-.0259 

-.0378 

-.0538 

it. 

.0101 


.06831 

.3500 

-.0683 

-.0252 

-.0101 

-.0468 

-.0420 

-.0558 

— 025 Q 

-.0250 

-.0558 

-.0420 

5 

-.0054 

-.0101 

-.0232 

-.0683 

.2500 

.0683 

.0232 

-.0558 

-.0378 

-.0259 

-.0180 

-.0293 

-.0339 

-.0359 

6 

-.0156 

-.0218 

-.0378 

-.0683 

.0683 

.2500 

.0683 

-.0683 

-.0378 

-.0218 

-.0136 

-.0339 

-.0359 

-.0339 

7 

-.0180 

-.0259 

-.0378 

-.0538 

.0232 

,0683 

.2500 

-.0683 

-.0232 

-.0101 

-.0054 

-.0359 

-.0339 

-.0293 

8 

-.0259 

-.0358 

-.0420 

-.0467 

-.0101 

-.0252 

-.0683 

.2500 

.0683 

.0232 

.0101 

-.0420 

-.0338 

-.0259 

9 

-.0558 

-.0420 

-.0467 

-.0420 

-.0218 

-.0378 

-.0683 

.0683 

.2500 

.0683 

.0232 

-.0538 

-.0378 

-.0259 

10 

-.0420 

-.0467 

-.0420 

-.0338 

-.0259 

-.0378 

-.0538 

.0232 

.0683 

.2500 

.0683 

-.0683 

-.0378 

-.0218 

11 

-.0467 

-.0420 

-.0338 

-0259 

-.0259 

-.0338 

-.0420 

.0101 

.0232 

.0683 

.2500 

-.0683 

-.0232 

-.0101 

12 

-.0538 

-.0378 

-.0259 

-.0180 

-.0293 

-.0339 

-.0359 

-.0054 

-.0101 

-.0232 

-.0683 

.2500 

.0683 

.0232 

15 

-.0685 

-.0578 

-.0218 

-.0136 

-0339 

-.0359 

-.0339 

-.0136 

-.0218 

-.0378 

-.0683 

.0683 

.2500 

.0683 

14 

-.0683 

-.0232 

-.0101 

-.0054 

-.0359 

-0339 

-.0293 

-.0180 

-.0259 

-.0378 

-.0538 

.0232 

.0683 

.2500 


•p- 

VO 


KACA TN h086 




50 


NACA TN 14086 


TABLE HI.- APPROXIMATE VALUES OF M FOR DIRICHLET PROBLEM 


M 

Iteration number 

0 

B 

2 

3 

B 

5 

6 

B 

8 

9 

10 

1 

1.270 

1.688 

1.767 

1.762 

1.743 

1.727 

1.717 

1.711 

1.708 

1.706 

1.705 

2 

2.170 

2.586 

2.613 

2.579 

2.547 

2.526 

2.513 

2.505 

2.501 

2.499 

2.497 

5 

-1.970 

-2.541 

-2.775 

-2.884 

-2.939 

-2.967 

-2.982 

-2.991 

-2.995 

-2.997 

-2.999 

4 

-1.520 

-2.182 

-2.469 

-2.592 

-2.647 

-2.674 

-2.687 

-2.694 

-2.698 

-2.770 

-2.701 

5 

3.360 

4.269 

4.508 

4.574 

4.592 

4.596 

4.598 

4.598 

4.598 

4.597 

4.597 

6 

.84o 

I.O83 

1.156 

1.183 

1.193 

1.198 

1.200 

1.200 

1.201 

1.201 

1.201 

7 

-It-. 100 

-4.925 

-5.086 

-5.109 

-5.108 

-5.104 

-5.101 

-5.100 

-5.099 

-5.098 

-5.098 

8 

1.170 

1.961 

2.344 

2.525 

2.611 

2.653 

2.673 

2.684 

2.689 

2.692 

2.694 

9 

1.580 

2.409 

2.805 

2.999 

3.096 

3.145 

3.171 

3.184 

3.190 

3.194 

3.196 

10 

.900 

1.393 

1.663 

1.815 

1.899 

1.944 

1.968 

1.980 

1.987 

1.991 

1.995 

11 

-1.220 

-1.292 

-1.204 

-1.122 

-1.067 

-1.036 

-1.018 

-1.008 

-1.003 

-1.001 

-.999 

12 

-1.280 

-1.669 

-1.844 

-1.927 

-1.967 

-1.986 

-1.995 

-2.000 

-2.001 

-2.002 

-2.003 

13 

-,2ltO 

-.621 

-.865 

-991 

-1.052 

-1.081 

-1.095 

-1.101 

-1.104 

-1.106 

-1.106 

14 

-1.790 

-2.452 

-2.748 

-2.884 

-2.946 

-2.975 

-2.989 

-2.995 

-2.998 

-3.000 

-3.001 


TABLE IV.- ERROR OF APPROXIMATE SOLUTION FOR DIRICHLET PROBLEM 61 



Iteration number 

€ 

B 

2 . . 

3 

B 

5 

6 

7 

8 

9 

10 

1 

0.418 

0.080 

-0.006 

-0.019 

-0.015 

-0.010 

-0.006 

-0.003 

-0.002 

-0.001 

2 

.416 

.028 

-.034 

-.032 

-.021 

■ 

-.007 

-.004 

-.002 

-.001 

3 

-.571 

-.233 

-.109 

-.055 

-.029 

BKjuJ 

-.008 

-.004 

-.002 

-.001 

4 

-.662 

-.286 

-.123 

-.056 

-.027 


-.007 

-.004 

-.002 

-.001 

5 

.909 

.239 

.065 

.018 

.005 

.001 

0 

0 

0 

0 

6 

.243 

.073 

.027 

.011 

.004 

.002 

.001 

0 

0 

0 

7 

-.825 

-.161 

-.023 

.002 

.004 

.003 

.001 

.001 

0 

0 

8 

.791 

• 383 

.180 

.086 

.042 

.021 

• Oil 

.005 

• 

O 

O 

.001 

9 

.829 

.396 

.194 

.097 

.049 

.025 

.013 

.007 

.004 

.002 


.493 

.270 

.152 

.084 

.045 

.024 


.007 

.004 

.002 


-.072 

.088 

.082 

.054 

.032 

.018 

.010 

.005 

.003 

.001 

12 

-.389 

-.175 

-.083 

-.040 

-.019 

-.009 

-.oo4 

-.002 

-.001 

0 

13 

-.381 

-.244 

-.126 

-.061 

-.029 

-.014 

-.006 

-.003 

-.001 

-.001 

14 

-.662 

-.297 

-.135 

-.062 

-.029 

-.014 

-.006 

-.003 

-.001 

-.001 


a £ (k+l) 


is error of kfch approximation 
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TABLE V.- MATRICES FOR NEUMANN PROBLEM - Continued 


0 3 ) (W) t ^rix 
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TABLZ V.- MATRICES FOB NEUMAHN FHOBIS4 - Continued 


(a) A r matrix 
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TABIE VI.- APPROXIMATE VALUES CF m FOR KEIMASH FROBIH1 

[oiven values of Ajp taken as are -0.60, -0.1(6, 0.32, -1.06, 1.48, 

0.65, -0.03, -0.64, -0.42, -0.77, -O.78, 1.06, 0.30, and O.95] 




Iteration number 


1 

2 

5 

4 

5 

6 

7 

8 

9 

10 

11 

12 

15 

1* 

1 

-0.956 

-1.163 

-1.278 

-1.338 

-1.365 

-1.367 

-1.359 

-1.341 

-1.319 

-1.293 

-I.265 

-1.235 

-1.2C5 

-1.175 

2 

-767 

-979 

-1.129 

-1.241 

-1.327 

-1.395 

-1.449 

-1.493 

-1.529 

-1.558 

-1.581 

-1.599 

-1.613 

-1.623 

3 

.674 

1.020 

1.338 

1.621 

1.866 

2.077 

2.257 

■'Iff 1 

2.539 

2.649 

2.743 

2.824 

2.892 

2.952 

4 

-1.178 

-2.286 

-2.647 

-2.910 

-3-104 

-3.247 

-3-353 

R 1 ill 

-3.487 

-3.526 

-3.552 

-3.567 

-3.574 

-5-575 

5 

2.466 

3.148 

3.637 

3.999 

4.273 

4.485 

4.653 

W 1 !K1 

4.898 

4.989 

5.066 

5.132 

5.189 

5-24o 

6 

.960 

1.152 

1.246 

1.301 

1.338 

1.366 

1.390 

1.4-12 

1.432 

1.452 

1.471 

1.489 

1.507 

1.525 

7 

-.235 

-.473 

-.696 

-.881 

-1.027 

-1.136 


ESB3 

-1.305 

-1-325 

-1.554 

-1.334 

-1.328 

-1.517 

8 

-.964 

-1.179 

-1-297 

-1.374 

-1.426 

-1.461 

-1.486 

-1.501 

-1.509 

-1.511 

-1.506 

-1.501 

-1.491 

-1.477 

9 

-.568 

--587 

-.548 

-.486 

-.416 

-.348 

-.283 

-.224 

-.170 

-.122 

-.078 

-.038 


-.031 

10 

-1.195 

-1.435 

-1.573 

-1.655 

-1.706 

-1.737 

-1.758 

-1.771 

-1-779 

-1.783 

-1.785 

-1.785 

-1.785 

-1.780 

11 

-1.181 

-1.373 

-1.451 

-1.466 

-1.448 

-1.414 

-1.371 

-1-327 

-1.262 

-1.238 

-I.196 

-1.156 

-1.118 

-1.082 

12 

1.863 

2.478 

2.955 

3.327 

3.621 

3-855 

4.043 

4.194 

4.317 

4.418 

4.500 

4.569 

4.626 

4.674 

13 

.320 

2.052 

.034 

-.149 

-325 

-.483 

-.619 

-733 

-.826 

-.904 

-965 

-1.013 

-1.050 

-1.077 

14 

1.537 

1.920 

2.178 

2.364 

2.507 

2.625 

2.722 

2.808 

2.886 

2.956 

3.021 

3.061 

3.136 

3.188 


TABIB VII— ERROR OP AEEROmAIB SQUHIOH FOR KEUMAKK PROBLEM 


■ 

Iteration number 


1 

2 

3 

4 

5 

6 


8 

9 

IO 

11 

12 

13 

14 

H 

-0.356 

-0.207 

-0.115 

-0.060 

-0.025 

-0.004 

0.009 

0.017 

0.023 

0.026 

0.028 

0.029 

0.030 

0.030 

2 

-.307 

-.211 

-.151 

-.112 

-.086 

-.068 

-.054 

-.044 

-.056 

-.269 

-.231 

-.018 

—oi4 

-.010 

3 

.354 

.346 

.318 

.282 

.245 

.211 

.280 

• 153 

.230 

.110 

.094 

.080 

.069 

.059 

4 

—722 

-.505 

-.361 

-.263 

-.194 

-.143 

-.106 

-.078 

-.056 

-.039 

— 026 

-.015 

-.007 

-.001 

5 

.gee 

.682 

.489 

.362 

.274 

.212 

.168 

.135 

.110 

.091 

.077 

.066 

.057 

.051 

6 

.330 

.172 

.094 

.055 

.037 

.028 

.024 

.022 

.020 

.020 

.019 

.018 

.018 

.017 

7 

-.203 

-.241 

-.223 

-.185 

-.145 

-.109 

-.079 

-055 

-.035 

-.020 

-.009 

0 

.006 

.011. 

8 

-.344 

-.195 

-.118 

-.077 

-.052 

.036 

-.024 

-.015 

-.008 

-.002 

.003 

.007 

.011 

.014 

9 

-.148 

-.019 

.039 

.063 

.069 

.069 

.064 

.059 

.054 

.049 

.044 

.040 

.036 

.053 

10 

-.425 

—240 

-.138 

-.082 

-.050 

-.052 

-.020 

-.013 

-.008 


-.002 

0 

.002 

.005 

11 

-.401 

-.192 

-077 

-.015 

.018 

.035 

.042 

.045 

.045 

.044 

.042 

.o4o 

.038 

.036 

12 

.803 

.615 

.476 

.372 

.294 

.234 

.188 

.151 

.123 

.101 

.085 

.068 

-057 

.048 

13 

.020 

-.115 

-.171 

-.184 

-.176 

-.158 

-.156 

-.115 

-.094 

-077 

— 06l 

-.048 

—036 

—027 

1* 

.590 

• 380 

.258 

.186 

.143 

.ns 

.099 

.087 

.078 

.071 

.065 

.060 

.056 

.052 


a (k+l) 


la error of kth approximation. 
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MCA IN 4o8 6 



• Inner point 

A Normal point (l to 10 j v= 10) 

□ Boundary point (l to 14; r = 14) 

- - Normal segment (unit length) 

(1 to 14; 3 = 14) 

Figure 1.- Region shewing "boundary points, normal points, and inner 
normal segments used in calculation of C for Dirichlet problem. 
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□ Boundary point (l to 14 j y = 14) 

O Exterior normal point ( 15 to 32 ; £ = 32) 

— »■— Directed exterior normal segment (unit 
length) (l to 32; T = 32) 

Figure 2.- Region showing boundary points, exterior normal points, and 
exterior normal segments used in calculation of C* for Heumann 
problem. 
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— •» x-axis 


□ Boundary point (l to 14 5 y = 14) 

A Normal points (l to 10; V= 10) 

— Inner normal directed from normal to 
inner points (l to 14; p = 14) 


Figure 5.- Region showing boundary points, inner normal points, and inner 
segments used in calculation of C in Dirichlet problem. 
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Figure 5-- 

[eM] 



Number of approximation 


Maximum error 


e( k ) 


max 


max 


max 


1 ^ r £ Ik 


.(k+i>| 


at boundary for Neuman problem where 
and e r ( k+1 ) = AjU r - A^C 1 *) . 
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